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Summary 

A computer program called INHYD has been developed for intraply hybrid composite design. This 
report is a users manual for INHYD. INHYD embodies several composite micromechanics theories, 
intraply hybrid composite theories, and an integrated hygrothermomechanical theory. INHYD can 
be run in both interactive and batch modes. It has considerable flexibility and capability, which the 
user can exercise through several options. These options are demonstrated through appropriate 
INHYD runs in the manual. 


Introduction 

Intraply hybrid composites have been investigated theoretically and experimentally at Lewis over 
the past several years. Theories developed during these investigations and corroborated by attendant 
experiments were used to develop a computer program identified as INHYD (INtraply HYbrid 
composite Design). INHYD was designed to operate in the interactive and batch modes. This report 
is a users manual for INHYD. 

INHYD includes several composite micromechanics theories, intraply hybrid composite theories, 
and an integrated hygrothermomechanical theory. Equations from these theories are used by the 
program as appropriate for the user’s specific applications. The elements and subroutines in INHYD 
are summarized in appendix A. The capabilities of the INHYD computer code are described in 
appendix B. The symbols used in the program are defined in appendix C. Before running INHYD for 
the first time, the user should read the program capability description in appendix B, which is a 
reproduction of NASA TM-82593 included herein for completeness. 

First, the manual presents examples for the interactive mode — for a dry, room-temperature case 
and then for a case with moisture present. Then an interactive case is described, where the data are 
preentered as they would be for a batch case. Next batch cases are described. After following through 
the test cases given, a user should be able to use INHYD either in the interactive or batch mode and 
with or without temperature and moisture effects to predict properties for uniaxial intraply hybrid 
composites. 


INHYD Interactive Mode 

Ambient Conditions, Dry 

To run INHYD, the user must first install and compile the program on his/her computer according 
to the system in use. The procedure used on the Lewis Research Center’s UNIVAC 1 100 is described 
in detail here by using the actual computer printout. Computer prompt signals are identified with 
upper-case letters. User entries follow the computer symbol >. Appropriate explanatory notes 
follow each printout. 


>(?3sar3 inhyd. 

READY 

>0ni3Pfi *3bs 

HAP 30RKC1) S74T11 04/22/83 11H2:46 (->0) 
>in inhyd* 


>end 

START=01251Af PROG SI ZE ( I /D > =1 0822/4956 
SYS$*RLIBi. LEVEL 
END MAP* errors: 0 
>0KOt *3bs 


INHYD TO BE RUN BATCH OR INTERACTIVE <B/I) 



As illustrated, the user is prompted by the computer, and his/her responses or entries are the input 
required for the run. 

INHYD INTERACTIVE DATA ENTRY 

SECONDARY COMPOSITE VOLUME RATIO 

> ♦ 2 

The user selected 0.2 for the secondary composite volume ratio. 

WHICH COMPOSITE DATA TO BE ENTERED (P/S) 

>p 

The choice was primary (P) or secondary (S); the user selected primary. 

PROPERTIES TO BE ENTERED (C/F/M/T) 

>f 

The choice was composite (C), fiber (F), dry matrix (M), or matrix with moisture and temperature 
effects (T); the user selected fiber properties for the primary composite. 

INPUT DATA HEADING ( 1 TO 80 CHARACTERS) 

> 3 s Graphite 

This heading identifies the primary fiber property table to be printed out by the computer. 


ENTER THE FIBER PROPERTY VALUES 

EFPl EFP2 GFP12 GFP23 NUFP12 NUFP23 

>*32e8T *2e7» .237, . Ie7f *2» *25 

CTEFPl CTEFP2 RHOFP NFP DIFP CFPC 

>“.56e“6T «56e — 5f *63e“lf *le5r *3e — 3r *17 

KFPl KFP2 KFP3 SFPT SFPC 

>*58e3» ♦58e2r ♦58e2> *Ae6f *4e6 

Here the user has entered the 17 values for the primary fiber properties that will be used in the 
micromechanics for calculating the primary composite properties. (Entry headings are defined in 
appendix C.) 

PROPERTIES TO BE ENTERED (C/F/M/T) 

>m 

Dry matrix (M) was selected. 

INPUT DATA HEADING ( 1 TO 80 CHARACTERS) 

>3501 -5 T dry f 70 decrees f 

This heading identifies the primary matrix property table to be printed out by the computer. 


ENTER THE MATRIX PROPERTY VALUES 
EMP GMP NUMP CTEMP RHOMP CMPC 

>.4Ae6f * 1643e6» .A? *32e-4» *44 3e-l . .25 

KMP SMPT SMPC SMPS BTAMP DIFMP 

>.125el f .68e4f .363e5f .7e4T .4, *2e-3 

Here the user has entered the 12 values for the primary matrix properties that will be used in the 
micromechanics for calculating the primary composite properties. (Entry headings are defined in 
appendix C.) 


COMPOSITE HEADING ( 1 TO 80 CHARACTERS) 



This heading identifies the primeiry composite property table. Table entries are to be calculated from 
the primeiry fiber properties and primary matrix properties and tabulated under this heading. The 
primary composite is AS/E (graphite/epoxy). 

ENTER CONSTITUENT VOLUME RATIOS 
FIBER MATRIX VOID 
>*55r *AAf .01 

The fiber volume ratio, matrix volume ratio, and void volume ratio of the primary composite are 
entered as desired by the user. The sum of the three values must be 1.0. 

VOID CONDUCTIVITY 
> ♦ 1 

The thermal conductivity of the voids was chosen as 0.1. 

METHOD FOR CALCULATING COMPOSITE PROPERTIES <H/M) 

>m 

The user can choose hygrothermomechanical theory (H) or micromechanics theory (M) for 
calculating the primary composite properties. Micromechanics theory was selected since moisture 
and temperature effects are not being studied. 

WHICH COMPOSITE DATA TO BE ENTERED (P/S) 

>s 

Now the secondary composite (S) properties will be calculated from the fiber and matrix properties 
about to be entered below. Since the procedure is the same as that used for the primary composite, 
fewer comments will be required. 

PROPERTIES TO BE ENTERED (C/F/M/T) 

>f 


INPUT DATA HEADING ( 1 TO 80 CHARACTERS) 

> S - d 1 3 S S 

The secondary fiber is S-Glass. 


ENTER THE FIBER PROPERTY VALUES 

EFSl EFS2 GFS12 GFS23 NUFS12 NUFS23 

>.12AeQf .124e3f .517e7T .517e7» .2 f .2 

CTEFSl CTEFS2 RHOFS NFS DIFS CFSC 

>.28e-Sf .28e-5, .9e-lf .204e3F .36e-3r .17 

KFSl KFS2 KFS3 SFST SFSC 

>.75el r *75el r .75elf * 36e 6 r *3e6 


PROPERTIES TO BE ENTERED (C/F/M/T) 

> III 


INPUT DATA HEADING ( 1 TO 80 CHARACTERS) 
>3501 ~5 T dry f 70 decrees f 


ENTER THE MATRIX PROPERTY VALUES 
EMS GMS NUMS CTEMS RHOMS CMSC 

>*A6g6, .1643e6f .4 f . 32e-4f .443e-l t .25 

KMS SMST SMSC SMSS BTAMS DIFMS 

>.125elF *68eAf .363e5F *7e4f .4 f . 2e-3 

The same matrix is used for the secondary composite as was used for the primary one. If desired, a 
different matrix could be used. 



COMPOSITE HEADING (1 TO 80 CHARACTERS) 

>5-dl3ss/e 


ENTER CONSTITUENT VOLUME RATIOS 
FIBER MATRIX VOID 
>*55f ♦44f .01 


VOID CONDUCTIVITY 

>.l 


METHOD FOR CALCULATING COMPOSITE PROPERTIES (H/M) 


HYBRID HEADING ( 1 TO 80 CHARACTERS) 

>3s/e//s-i3lass/ey 80/20 

This is the heading for the hybrid composite table that will be calculated. 

METHOD FOR CALCULATING HYBRID PROPERTIES (R/M) 

>r 

The choice is micromechanics theory (M) or rule of mixtures (R); the user selected the rule of 
mixtures. 

The computer then lists the headings and properties previously entered by the user for the fibers (17 
properties) and matrix (12 properties) of the primary composite in tabular form. Next it lists the 
respective volume ratios and the 37 calculated properties for the primary composite and a table of 
flexural strengths for the primary composite as calculated by several theories. 


PRIMARY FIBER PROPERTIES? AS GRAPHITE 


1 

ELASTIC MODULI 

EFPl 

.3200+08 

2 


EFP2 

.2000+07 

3 

SHEAR MODULI 

GFP12 

.2000+07 

A 


GFP23 

.1000+07 

5 

POISSON'S RATIO 

NUFP12 

.2000+00 

6 


NUFP23 

.2500+00 

7 

THERM. EXP. COEF. 

CTEFPl 

-.5600-06 

8 


CTEFP2 

.5600-05 

9 

DENSITY 

RHOFP 

♦6300-01 

10 

NO. OF FIBERS/END 

NFP 

.1000+05 

11 

FIBER DIAMETER 

DIFP 

.3000-03 

12 

HEAT CAPACITY 

CFPC 

.1700+00 

13 

HEAT CONDUCTIVITY 

KFPl 

♦5800+03 

14 


KFP2 

.5800+02 

15 


KFP3 

♦5800+02 

16 

STRENGTHS 

SFPT 

.4000+06 

17 


SFPC 

.4000+06 


RIMARY 

MATRIX PROPERTIES? 3501 

-5f DRY»70 

DEGREES F 

1 

ELASTIC MODULUS 

EMP 

.4600+06 

2 

SHEAR MODULUS 

GMP 

.1643+06 

3 

POISSON'S RATIO 

NUMP 

.4000+00 

4 

THERM. EXP. COEF. 

CTEMP 

.3200-04 

5 

DENSITY 

RHOMP 

.4430-01 

6 

HEAT CAPACITY 

CMPC 

.2500+00 

7 

HEAT CONDUCTIVITY 

KMP 

.1250+01 

8 

STRENGTHS 

SMPT 

.6800+04 

9 


SMPC 

.3630+05 

10 


SMPS 

.7000+04 

11 

MOISTURE COEF 

STAMP 

.4000+00 

12 

DIFFUSIVITY 

DIFMP 

.2000-03 


PRIMARY COMPOSITE PROPERTIES? AS/E 

BASED ON MIGROMECHANICS OF INTRAPLY HYBRID COMPOSITES: 
ELASTIC AND THERMAL PROPERTIES. 



FIBER VOLUHE RATIO - 
VOID CONDUCTIVITY - 


.550 

♦ 1000+00 


HATRIX VOLUHE RATIO - .440 

VOID VOLUME RATIO - .010 


1 

ELASTIC MODULI 

EPCl 

.1784+08 

2 


EPC2 

♦1187+07 

3 


EPC3 

.,1187 + 07 

4 

SHEAR MODULI 

6PC12 

.7087+06 

5 


GPC23 

.5612+06 

6 


GPC13 

.7087+06 

7 

POISSON'S RATIO 

NUPC12 

.2900+00 

8 


NUPC23 

.3175+00 

9 


NUPC13 

.2900+00 

10 

THERM. EXP. COEF. 

CTEPCl 

-.1280-06 

11 


CTEPC2 

.,1563-04 

12 


CTEPC3 

.1563-04 

13 

DENSITY 

RHOPC 

.5414-01 

14 

HEAT CAPACITY 

CPC 

.,1988 + 00 

15 

HEAT CONDUCTIVITY 

KPCl 

♦3195+03 

16 


KPC2 

♦3702+01 

17 


KPC3 

.3702+01 

18 

STRENGTHS 

SPCIT 

.2230+06 

19 


SPCIC 

.2230+06 

20 


SPC2T 

.4944+04 

21 


SPC2C 

.2639+05 

22 


SPC12 

.5335+04 

23 

MOIST. DIFFUSIVITY 

DPCl 

.8800-04 

24 


DPC2 

.5168-04 

25 


DPC3 

.5168-04 

26 

MOIST. EXP. COEF. 

BTAPCl 

.5402-02 

27 


BTAPC2 

♦1441+00 

28 


BTAPC3 

.1441+00 

29 

FLEXURAL MODULI 

EPCIF 

. 1784 + 08 

30 


EPC2F 

.,1187 + 07 

31 

STRENGTHS 

SPC23 

.5868+04 

32 


SPCIF 

.2788+06 

33 


SPC2F 

. 1041 + 05 

34 


SPCSB 

.8002+04 

35 

PLY THICKNESS 

TPC 

.5000-02 

36 

INTERPLY THICKNESS 

PLPC 

.,5850-04 

37 

INTERFIBER SPACING 

PLPCS 

.,5850-04 




FLEXURAL STRENGTH 


I THEORY I LONGITUDINAL I TRANSVERSE I 

I BILINEAR I .2230+06 I .8328+04 I 

I SOLUTION III 

I LINEAR RECTANGULAR I .2230+06 I .9737+04 I 

I SOLUTION III 

I NEUTRAL AXIS SHIFT I .2230+06 I .9737+04 I 

I SOLUTION III 

I PARABOLIC I .2788+06 I .1041+05 I 

I SOLUTION III 

I RECTANGULAR I .3345+06 I .1249+05 I 

I SOLUTION I I I 


The longitudinal and transverse flexural strengths listed in the composite property table (items 32 and 
33) are calculated by the parabolic theory. The ply thickness (item 35) is internally set at 0.005 in., 
which is a typical value for most fiber composites. The ply thickness is used to calculate flexural 
strength. 

The computer repeats this procedure for the secondary composite — first listing the data entered by 
the user and then tabulating the 37 calculated properties for the secondary composite (which in this 
example is S-Glass/epoxy) and the flexural strengths. 



SECONDARY FIBER PROPERTIES? S-GLASS 


1 

ELASTIC MODULI 

EFSl 

.1240+08 

2 


EFS2 

*1240+08 

3 

SHEAR MODULI 

GFS12 

*5170+07 

4 


GFS23 

♦5170+07 

5 

POISSON'S RATIO 

NUFS12 

*2000+00 

6 


NUFS23 

*2000+00 

7 

THERM* EXP* COEF* 

CTEFSl 

*2800-05 

8 


CTEFS2 

*2800-05 

9 

DENSITY 

RHOFS 

*9000-01 

10 

NO* OF FIBERS/END 

NFS 

*2040+03 

11 

FIBER DIAMETER 

DIFS 

*3600-03 

12 

HEAT CAPACITY 

CFSC 

.1700+00 

13 

HEAT CONDUCTIVITY 

KFSl 

*7500+01 

14 


KFS2 

*7500+01 

15 


KFS3 

*7500+01 

16 

STRENGTHS 

SFST 

*3600+06 

17 


SFSC 

*3000+06 


SECONDARY 

MATRIX PROPERTIES? 3501 

-5f DRYf70 

DEGREES F 

1 

ELASTIC MODULUS 

EMS 

*4600+06 

2 

SHEAR MODULUS 

6MS 

.1643+06 

3 

POISSON'S RATIO 

NUMS 

.4000+00 

4 

THERM. EXP. COEF. 

CTEMS 

.3200-04 

5 

DENSITY 

RHOMS 

*4430-01 

6 

HEAT CAPACITY 

CMSC 

*2500+00 

7 

HEAT CONDUCTIVITY 

KMS 

.1250+01 

8 

STRENGTHS 

SMST 

*6800+04 

9 


SMSC 

.3630+05 

10 


SMSS 

.7000+04 

11 

MOISTURE COEF 

BTAMS 

.4000+00 

12 

DIFFUSIVITY 

DIFMS 

.2000-03 


SECONDARY 

COMPOSITE PROPERTIES? S-GLASS/E 


BASED 

ON MICROMECHANICS OF INTRAPLY HYBRID 

COMPOSITES 

ELASTIC AND THERMAL PROPERTIES 

• 


FIBER 

VOLUME RATIO - .550 

MATRIX VOLUME RATIO - 

VOID 

CONDUCTIVITY - *1000+00 

VOID VOLUME 

RATIO - . 

1 

ELASTIC MODULI 

ESCl 

.7061+07 

2 


ESC2 

*1881+07 

3 


ESC3 

.1881+07 

4 

SHEAR MODULI 

6SC12 

*8449+06 

5 


GSC23 

*8449+06 

6 


6SC13 

,8449+06 

7 

POISSON'S RATIO 

NUSC12 

.2900+00 

8 


NUSC23 

.2900+00 

9 


NUSC13 

.2900+00 

10 

THERM. EXP. COEF. 

CTESCl 

.3779-05 

11 


CTESC2 

.1348-04 

12 


CTESC3 

.1348-04 

13 

DENSITY 

RHOSC 

.6899-01 

14 

HEAT CAPACITY 

CSC 

.1926+00 

15 

HEAT CONDUCTIVITY 

KSCl 

.4675+01 

16 


KSC2 

.2750+01 

17 


KSC3 

.2750+01 

18 

STRENGTHS 

SSCIT 

.2050+06 

19 


SSCIC 

. 1708+06 

20 


SSC2T 

.4223+04 

21 


SSC2C 

.2254+05 

22 


SSC12 

♦4842+04 

23 

MOIST. DIFFUSIVITY 

DSCl 

*8800-04 

24 


DSC2 

*5168-04 

25 


DSC3 

*5168-04 

26 

MOIST* EXP* COEF* 

BTASCl 

*1365-01 

27 


BTASC2 

*1433+00 

28 


BTASC3 

*1433+00 

29 

FLEXURAL MODULI 

ESCIF 

*7061+07 

30 


ESC2F 

*1881+07 


*440 

010 



31 

STRENGTHS 

SSC23 

.4842+04 

32 


SSCIF 

.2329+06 

33 


SSC2F 

.8892+04 

34 


SSCSB 

♦7263+04 

35 

PLY THICKNESS 

TSC 

.5000-02 

36 

INTERPLY THICKNESS 

PLSC 

^7020“04 

37 

INTERFIBER SPACING 

PLSCS 

♦7020-04 




I THEORY 


BILINEAR 

SOLUTION 


FLEXURAL STRENGTH 


LONGITUDINAL 


TRANSVERSE 


.1864+06 


.7114+04 


LINEAR RECTANGULAR 
SOLUTION 


.2019+06 


.8317+04 


1 «« :|c }|u|c« )|n|c ««««««« ;tnK «« )|c « )|c «« :|n|n|c 3|c « ]|n| }tc «« ]|nK «« 3|c ««« }|c «««««««« iK«« «« 1 






I NEUTRAL AXIS SHIFT I .2019+06 I .8317+04 I 

I SOLUTION III 

I PARABOLIC I .2329+06 I .8892+04 I 

I SOLUTION III 

I RECTANGULAR I .2795+06 I .1067+05 I 

I SOLUTION III 


Finally the calculated hybrid composite property (35 properties) and hybrid composite flexural 
strength tables are printed out. 


HYBRID COMPOSITE PROPERTIES? AS/E//A S-GLASS/E > 80/20 


BASED ON PREDICTION OF PROPERTIES OF INTRAPLY HYBRID COMPOSITES 
(RULE OF MIXTURES). 


PRIMARY COMPOSITE VOLUME RATIO - .800 

SECONDARY COMPOSITE VOLUME RATIO - .200 


1 

ELASTIC MODULI 

EHCl 

♦1568+08 

2 


EHC2 

.1326+07 

3 


EHC3 

.1326+07 

4 

SHEAR MODULI 

6HC12 

.7360+06 

5 


GHC23 

.6180+06 

6 


GHC13 

.7360+06 

7 

POISSON'S RATIO 

NUHC12 

.2900+00 

8 


NUHC23 

.3120+00 

9 


NUHC13 

.2900+00 

10 

THERM. EXP. COEF. 

CTEHCl 

♦6534-06 

11 


CTEHC2 

♦1520-04 

12 


CTEHC3 

.1520-04 

13 

DENSITY 

RHOHC 

.5711-01 

14 

HEAT CAPACITY 

CHC 

.1976+00 

15 

HEAT CONDUCTIVITY 

KHCl 

.2566+03 

16 


KHC2 

.3511+01 

17 


KHC3 

.3511+01 

18 

STRENGTHS 

SHCIT 

.2194+06 

19 


SHClC 

.2126+06 

20 


SHC2T 

.4800+04 

21 


SHC2C 

.2562+05 

22 


SHC12 

.5236+04 

23 

MOIST. DIFFUSIVITY 

DHCl 

.8800-04 

24 


DHC2 

.5168-04 

25 


DHC3 

.5168-04 

26 

MOIST. EXP. COEF. 

BTAHCl 

♦7052-02 

27 


BTAHC2 

♦1440+00 

28 


BTAHC3 

.1440+00 

29 

FLEXURAL MODULI 

EHCIF 

.1568+08 



30 


EHC2F 

.1326+07 

31 

STRENGTHS 

SHC23 

.5663+04 

32 


SHCIF 

.2699+06 

33 


SHC2F 

.1011+05 

34 


SHCSB 

.7854+04 

35 

FIBER VOL. RATIO 

VFH 

.5500+00 




I THEORY 


I BILINEAR 
I SOLUTION 


FLEXURAL STRENGTH 


LONGITUDINAL 


TRANSVERSE 


.2159+06 


.8085+04 


I LINEAR RECTANGULAR 
I SOLUTION 


.2193+06 


.9453+04 


NEUTRAL AXIS SHIFT 
SOLUTION 


.2193+06 


.9453+04 




I }|c 4c )K K « }K « )|c )|c 3K }|n(c}f( }|[ :(c :K }|c }|c :|c ]|c « }|nK « « t « )|nK ^ ^ }|( « ))c :tc 3K }|c ^ * )|c )K ^ ^ ^ ^ :1c ^ :K ^ * lc 4c « « « 3^ « I 




1 4c 4c 4c 4c 4c 4( 4c 4C 4c 4C 4c 4c 4( 4C «« 4c 4c 4c 4c 4c 4c 4c 4C 4c 4t 4c 4c 4c 4c 4( 4c 4( 4c 4c 4c 4c 4c 4c 4c 4c 4( 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4cl 


I PARABOLIC I .2699+06 I .1011+05 I 

I SOLUTION III 

1 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c ic 4c 4c 4c 4c 4c 4c 4c 4c 4c I 
I RECTANGULAR I ,3239+06 I ,1213+05 I 

I SOLUTION I I I 


With Moisture Effects 


An example for the same composite, 80/20 AS/E//S-G/E, but with moisture and temperature 
effects is presented below. Only features differing from those explained in the previous “dry” 
example are explained. 


>0asdfd inhvfd. 

READY 

>0m3Pri .abs 

MAP 30RKC1) S74T11 04/22/83 11:37:55 (->0) 
>in inhyd. 

>end 

START=012516f PROG SIZE( I/D)=10822/4956 
SYS$4cRLIB$. LEVEL 
END MAP. errors: 0 
>0xat .abs 

INHYD TO BE RUN BATCH OR INTERACTIVE (B/I) 
>i 


INHYD INTERACTIVE DATA ENTRY 

SECONDARY COMPOSITE VOLUME RATIO 

>.2 

WHICH COMPOSITE DATA TO BE ENTERED (P/S) 

>p 


PROPERTIES TO BE ENTERED (C/F/M/T) 
>f 


INPUT DATA HEADING < 1 TO 80 CHARACTERS) 
>as draphite 



ENTER THE FIBER PROPERTY VALUES 
EFPl EFP2 GFP12 GFP23 NUFP12 NUFP23 

>*32e8f *2e7f *2e7j *le7T •2> *25 

CTEFPl CTEFP2 RHOFP NFP DIFP CFPC 

>-*56e-Af .56e-5» .63e-lf le5r *3e“3> *17 

KFPl KFP2 KFP3 SFPT SFPC 

> « 58e3 r *5862 > ♦ 58e2 f ♦ 4e6 f ♦ Ae6 


PROPERTIES TO BE ENTERED (C/F/M/T) 

>t 

The matrix with moisture and temperature properties (T) was selected. Now additional input is 
needed for the matrix. Nineteen entries are needed instead of the 12 required for the previous case: 

INPUT DATA HEADING ( 1 TO 80 CHARACTERS) 

>3501-5fwith 2% moisture 


ENTER THE MATRIX PROPERTY VALUES 
EMP GMP NUMP CTEMP RHOMP CMPC 

>.46eAf .lA43e6» .4 t . 32e-4f .443e-l r .25 

KMP SMPT SMPC SMPS STAMP DIFMP 

>*125el 9 *68e4f «363e5? .7e4» *4> *2e-3 


ENTER THE THERMAL/MOISTURE PROPERTY VALUES 
TO TGDR TEMP MOIST (%) 

>273. f 460. f 294. f2.0 

The first three new entries in the above line are the reference temperature, the glass transition 
temperature of the resin, and the use temperature. Here they are entered in Kelvin, but any scale can 
be used since it is a ratio of differences. If Fahrenheit were used, the entries would be 32®, 368®, and 
70®, respectively. The last entry is for 2 percent moisture. In addition the following three entries are 
needed: 

KMST CMST RHOMST 

>.lf .1>1.0 

These entries are the conductivity, the heat capacity, and the density of the moisture substance (see 
INHYD.INSTR). 

COMPOSITE HEADING (1 TO 80 CHARACTERS) 

>3s/e with 2% moisture in matrix 


ENTER CONSTITUENT VOLUME RATIOS 
FIBER MATRIX VOID 
> * 55 f *AAf .01 


VOID CONDUCTIVITY 
>♦1 


METHOD FOR CALCULATING COMPOSITE PROPERTIES <H/M) 

>h 

Hygrothermomechanical theory was selected since 2 percent moisture is present in the matrix. 


UHICH COMPOSITE DATA TO BE ENTERED (P/S) 

>s 


PROPERTIES TO BE ENTERED (C/F/M/T) 

>f 


INPUT DATA HEADING < 1 TO 80 CHARACTERS) 
>5~alsSS 


ENTER THE FIBER PROPERTY VALUES 
EFSl EFS2 GFS12 GFS23 NUFS12 NUFS23 

>.124e8T *12468 f .517e7» ♦SI?©?. *2f *2 

CTEFSl CTEFS2 RHOFS NFS DIFS CFSC 

> ♦ 28e~5 . • 28 e- 5 f'*9e“lf * 204e3 . *366-3. *17 

KFSl KFS2 KFS3 SFST SFSC 

->*75el. *75el. *75el. *36e6. *3e6 


PROPERTIES TO BE ENTERED (C/F/M/T) 
>t 


INPUT DATA HEADING ( 1 TO 80 CHARACTERS) 
>3501-5 With 2% noisture 


ENTER THE MATRIX PROPERTY VALUES 
EMS GMS NUMS CTEMS RHOMS CMSC 

>*46e6. *1643e6. *4. *32e-4. *443e-l . *25 

KMS SMST SMSC SMSS BTAMS DIFMS 

>*125el. *68e4. «363e5. *7e4. *4. *2e— 3 


ENTER THE THERMAL/MOISTURE PROPERTY VALUES 
TO TGDR TEMP MOIST (%) 

>273* .460* ,294* .2.0 

KMST CMST RHOMST 

> * 1 . ♦ 1 . 1 * 0 


COMPOSITE HEADING ( 1 TO 80 CHARACTERS) 
>s-3l35s With 2 % moisture in mstriw 


ENTER CONSTITUENT VOLUME RATIOS 
FIBER MATRIX VOID 
>*55. *44. *01 


VOID CONDUCTIVITY 
x*l 


METHOD FOR CALCULATING COMPOSITE PROPERTIES (H/M) 
>h 


HYBRID HEADING ( 1 TO 80 CHARACTERS) 
>3s/e//s-d/e . 80/20 with 2X moisture in metrix 


METHOD FOR CALCULATING HYBRID PROPERTIES (R/M) 

>r 

The choice is rule of mixtures (R) or micromechanical theory (M); the user selected the rule of 
mixtures. 


PRIMARY FIBER PROPERTIES; AS GRAPHITE 
1 ELASTIC MODULI 


2 

3 SHEAR MODULI 

4 

5 POISSON'S RATIO 

6 

7 THERM* EXP* COEF* 


EFPl *3200+08 
EFP2 *2000+07 
GFP12 *2000+07 
GFP23 *1000+07 
NUFP12 *2000+00 
NUFP23 *2500+00 
CTEFPl -.5600-06 


10 



8 


CTEFP2 

.5600-05 

9 

DENSITY 

RHOFP 

.6300-01 

10 

NO. OF FIBERS/END 

NFP 

.1000+05 

11 

FIBER DIAMETER 

DIFP 

.3000-03 

12 

HEAT CAPACITY 

CFPC 

•1700+00 

13 

HEAT CONDUCTIVITY 

KFPl 

.5800+03 

14 


KFP2 

.5800+02 

15 


KFP3 

.5800+02 

16 

STRENGTHS 

SFPT 

.4000+06 

17 


SFPC 

.4000+06 

PRIMARY 

MATRIX PROPERTIES? 3501 

-SfWITH 2% 

MOISTURE 


REFERENCE TEHP. - 273*00 

TEST TEHP. - 294.00 

MOISTURE CONDUCTIVITY - .1000+00 

MOISTURE DENSITY - .1000+01 

DRY GLASS TRANS. TEMP. - 460.00 

PCT. MOISTURE - 2.000 

MOISTURE HEAT CAPACITY - .1000+00 


ORIGINAL MATRIX PROPERTIES 


1 

ELASTIC MODULUS 

EMP 

.4600+06 

2 

SHEAR MODULUS 

GMP 

. 1643+06 

3 

POISSON'S RATIO 

NUMP 

.4000+00 

4 

THERM. EXP. COEF 

. CTEMP 

.3200-04 

5 

DENSITY 

RHOMP 

.4430-01 

6 

HEAT CAPACITY 

CMPC 

.2500+00 

7 

HEAT CONDUCTIVITY KMP 

.1250+01 

8 

STRENGTHS 

SMPT 

.6800+04 

9 


SMPC 

.3630+05 

10 


SMPS 

♦7000+04 

11 

MOISTURE COEF 

BTAMP 

.4000+00 

12 

DIFFUSIVITY 

DIFMP 

.2000-03 

user’s input is tabulated in the preceding tables. 



NEW properties: 

CHANGED BECAUSE 



OF temperature/moisture effects 


t 

NEW EMP 

.3068+06 


t 

NEW GMP 

.1096+06 


t 

NEW NUMP 

♦4000+00 


t 

NEW CTEMP 

.4797-04 


* 

NEW RHOMP 

.4547-01 



NEW CMPC 

.3748+00 


» 

NEW KMP 

.1874+01 



NEW SMPT 

.4536+04 


t 

NEW SMPC 

.2421+05 



NEW SMPS 

.4669+04 


% 

NEW BTAMP 

. 4000+00 



NEW DIFMP 

.2000-03 



Temperature and moisture effects on the matrix are calculated and tabulated by the computer in the 
preceding table. These two tables are printed out side by side. 


PRIMARY COMPOSITE PROPERTIES? AS/E WITH 2X MOISTURE IN MATRIX 

BASED ON INTEGRATED THEORY FOR PREDICTING THE HYGROTHERMOMECHANICAL 

RESPONSE OF ADVANCED COMPOSITE STRUCTURAL COMPONENTS. 


FIBER 

VOLUME RATIO - 

.550 

MATRIX VOLUME RATIO 

VOID 

CONDUCTIVITY - ♦ 

1000+00 

VOID VOLUME 

RATIO - 

1 

ELASTIC MODULI 


EPCl 

♦1776+08 

2 



EPC2 

♦9275+06 

3 



EPC3 

♦9275+06 

4 

SHEAR MODULI 


6PC12 

.5181+06 

5 



6PC23 

.4346+06 



6 


GPC13 

.5181+06 

7 

POISSON'S RATIO 

NUPC12 

.2860+00 

8 


NUPC23 

.6711-01 

9 


NUPC13 

.2860+00 

10 

THERM. EXP. COEF. 

CTEPCl 

-.1208-06 

11 


CTEPC2 

.2146-04 

12 


CTEPC3 

.2146-04 

13 

DENSITY 

RHOPC 

♦5466-01 

14 

HEAT CAPACITY 

CPC 

♦2450+00 

15 

HEAT CONDUCTIVITY 

KPCl 

♦3198+03 

16 


KPC2 

.5407+01 

17 


KPC3 

.5407+01 

13 

STRENGTHS 

SPCIT 

♦2220+06 

19 


SPCIC 

.1231+06 

20 


SPC2T 

♦4470+04 

21 


SPC2C 

.2386+05 

22 


SPC12 

♦4601+04 

23 

MOIST. DIFFUSIVITY 

DPCl 

.8800-04 

24 


DPC2 

♦5168-04 

25 


DPC3 

♦5168-04 

26 

MOIST. EXP. COEF. 

BTAPCl 

.3620-02 

27 


BTAPC2 

♦ 1126 + 00 

28 


BTAPC3 

.1126+00 

29 

FLEXURAL MODULI 

EPCIF 

.1776+08 

30 


EPC2F 

.9275+06 

31 

STRENGTHS 

SPC23 

♦4601+04 

32 


SPCIF 

.1980+06 

33 


SPC2F 

.9411+04 

34 


SPCSB 

.6902+04 

35 

PLY THICKNESS 

TPC 

♦5000-02 

36 

INTERPLY THICKNESS 

PLPC 

.5850-04 

37 

INTERFIBER SPACING 

PLPCS 

.5850-04 




I THEORY 


I BILINEAR 
I SOLUTION 


FLEXURAL STRENGTH 


LONGITUDINAL 


TRANSVERSE 


♦1584+06 


.7529+04 






I LINEAR RECTANGULAR I .1937+06 I .8802+04 I 

I SOLUTION I I I 

1 NEUTRAL AXIS SHIFT I .1937+06 I .8802+04 I 

I SOLUTION I I I 

I PARABOLIC I .1980+06 I .9411+04 I 

I SOLUTION III 

I RECTANGULAR I .2376+06 I .1129+05 I 

I SOLUTION I I I 


The primary composite flexural strength table concludes the primary composite calculations. User 
input for the secondary composite is now tabulated. 


SECONDARY FIBER PROPERTIES? S-GLASS 


1 

ELASTIC MODULI 

EFSl 

. 1240+08 

2 


EFS2 

.1240+08 

3 

SHEAR MODULI 

GFS12 

.5170+07 

4 


GFS23 

.5170+07 

5 

POISSON'S RATIO 

NUFS12 

.2000+00 

6 


NUFS23 

♦2000+00 

7 

THERM. EXP. COEF. 

CTEFSl 

.2800-05 

8 


CTEFS2 

.2800-05 

9 

DENSITY 

RHOFS 

.9000-01 

10 

NO. OF FIBERS/END 

NFS 

♦2040+03 

11 

FIBER DIAMETER 

DIFS 

.3600-03 

12 

HEAT CAPACITY 

CFSC 

♦1700+00 

13 

HEAT CONDUCTIVITY 

KFSl 

♦7500+01 

14 


KFS2 

♦7500+01 

15 


KFS3 

♦7500+01 



STRENGTHS 


SFST 

SFSC 


3600+06 

3000+06 


16 

17 


SECONDARY MATRIX PROPERTIES? 3501-5 E WITH 2% MOISTURE 

REFERENCE TEMP. - 273.00 

TEST TEMP. “ 294.00 

MOISTURE CONDUCTIVITY - .1000+00 

MOISTURE DENSITY - .1000+01 

DRY GLASS TRANS. TEMP. - 460.00 

PCT. MOISTURE - 2.000 

MOISTURE HEAT CAPACITY - .1000+00 



ORIGINAL MATRIX PROPERTIES 


1 

ELASTIC MODULUS 

EMS 

♦4600+06 

2 

SHEAR MODULUS 

GMS 

.1643+06 

3 

POISSON'S RATIO 

NUMS 

.4000+00 

4 

THERM. EXPv COEF. 

STEMS 

.320C-04 

5 

DENSITY 

RHOMS 

. 4430-01 

6 

HEAT CAPACITY 

CMSC 

♦2500+00 

7 

HEAT CONDUCTIVITY 

KMS 

.1250+01 

8 

STRENGTHS 

SMST 

. 6800+04 

9 


SMSC 

.3630+05 

10 


3MSS 

.7000+04 

11 

MOISTURE COEF 

BTAMS 

.4000+00 

12 

DIFFUSIVITY 

DIFMS 

.2000-03 



NEW 

properties: 

CHANGED BECAUSE 


OF TEMPERATURE/MGISTURE EFFECTS 

* 

NEW 

EMS 

.3068+06 

* 

NEW 

GMS 

.1096+06 

* 

NEW 

NUMS 

. 4000+00 

* 

NEW 

CTEMS 

. 4797-04 

* 

NEW 

RHOMS 

♦ 4547-01 

* 

NEW 

CMSC 

.3748+00 

* 

NEW 

KMS 

.1874+01 


NEW 

oMST 

.4536+04 

* 

NEW 

SMSC 

. 2421+05 

* 

NEW 

3 MSS 

.4669+04 


NEW 

BTAMS 

♦ 4000 + 00 

* 

NEW 

D [FMS 

.2000-03 


SECONDARY COMPOSITE PROPERTIES? S-GLASS WITH 2% MOISTURE IN MATRIX 
BASED ON INTEGRATED THEORY FOR PREDICTING THE HYGROTHERMOMECHANICAL 
RESPONSE OF ADVANCED COMPOSITE STRUCTURAL COMPONENTS. 


FIBER VOLUME RATIO - .550 MATRIX VOLUME RATIO - .440 

VOID CONDUCTIVITY - .1000+00 VOID VOLUME RATIO - .010 


1 

ELASTIC MODULI 

ESC! 

.6981+07 

2 


ESC2 

.1303+07 

3 


ESC3 

. 1303+07 

4 

SHEAR MODULI 

GSC12 

.5872+06 

5 


GSC23 

.5872+06 

6 


GSC13 

.5872+06 

7 

POISSON'S RATIO 

NUSC12 

♦2860+00 

8 


NUSC23 

.1098+00 

9 


NUSC13 

.2860+00 

10 

THERM. EXP. COEF. 

CTESCl 

.3840-05 

11 


CTESC2 

.1932-04 

12 


CTESC3 

^1932-04 

13 

DENSITY 

RHOSC 

.6951-01 

14 

HEAT CAPACITY 

CSC 

.2289+00 

15 

HEAT CONDUCTIVITY 

KSCl 

.4950+01 

16 


KSC2 

.3617+01 

17 


KSC3 

.3617+01 


18 

STRENGTHS 

3SC1T 

.2027+06 

19 


SSCiC 

.9655+05 

20 


SSC2T 

.4470+04 

21 


SSC2C 

,2386+05 

22 


SSC12 

.4601+04 

23 

MOIST. DIFFUSIVITY 

DSCl 

.8800-04 

24 


DSC2 

.5168-04 

25 


DSC3 

.5168-04 

26 

MOIST. EXP. COEF. 

BTASCl 

.9209-02 

27 


BTASC2 

♦1102+00 

28 


BTASC3 

♦1102+00 

29 

FLEXURAL MODULI 

ESCIF 

.6981+07 

30 


ESC2F 

,1303+07 

31 

STRENGTHS 

SSC23 

.4601+04 

32 


SSCIF 

.1635+06 

33 


SSC2F 

.9411+04 

34 


SSCSB 

.6902+04 

35 

PLY THICKNESS 

TSC 

,5000-02 

36 

INTERPLY THICKNESS 

PLSC 

,7020-04 

37 

INTERFIBER SPACING 

PLSCS 

.7020-04 


I FLEXURAL STRENGTH I 
I *3(OK*»******#)K***)lc I 
I **#)!(*** ************** ********* ******#:<c#}|t#)<c********5*t **#*#*#*♦** I 
I THEORY I LONGITUDINAL I TRANSVERSE I 
!**************************************************************! 
I BILINEAR I .1308+06 I .7529+04 I 
I SOLUTION III 
I ******************************************* ***************)^*** I 
I LINEAR RECTANGULAR I .1650+06 I .8802+04 I 
I SOLUTION I I I 
!**************************************************************! 
I NEUTRAL AXIS SHIFT I .1650+06 I .8802+04 I 
I SOLUTION III 
!**************************************************************! 
I PARABOLIC I .1635+06 I .9411+04 I 
I SOLUTION I I I 
I ************************************************************** I 
I RECTANGULAR I *1962+06 I .1129+05 I 
I SOLUTION I I I 

I *********3MC** ************************************************* I 


HYBRID COMPOSITE PROPERTIES? AS/E//S-G/E t 80/20 WITH 2% MOISTURE IN MATRIX 
BASED ON PREDICTION OF PROPERTIES OF INTRAPLY HYBRID COMPOSITES 
(RULE OF MIXTURES/ ♦ 


PRIMARY COMPOSITE VOLUME RATIO - .800 

SECONDARY COMPOSITE VOLUME RATIO - »200 


1 

ELASTIC MODULI 

EHCl 

,1560+08 

2 


EHC2 

,1003+07 

3 


EHC3 

.1003+07 

4 

SHEAR MODULI 

GHC12 

.5319+06 

5 


GHC23 

.4651+06 

6 


GHC13 

.5319+06 

7 

POISSON'S RATIO 

NUHC12 

.2860+00 

8 


NUHC23 

.7566-01 

9 


NUHC13 

.2860+00 

10 

THERM. EXP. COEF, 

CTEHCl 

.6714-06 

11 


CTEHC2 

♦2103-04 

12 


CTEHC3 

.2103-04 

13 

DENSITY 

RHOHC 

.5763-01 

14 

HEAT CAPACITY 

CHC 

.2418+00 

15 

HEAT CONDUCTIVITY 

KHCl 

.2568+03 

16 


KHC2 

.5049+01 

17 


KHC3 

.5049+01 

18 

STRENGTHS 

SHCIT 

.2181+06 

19 


SHClC 

,1178+06 

20 


SHC2T 

.4470+04 


21 


SHC2C 

.2386+05 

22 


SHC12 

.4601+04 

23 

MOIST. DIFFUSIVITY 

DHCl 

.8800-04 

24 


DHC2 

.5168-04 

25 


DHC3 

.5168-04 

26 

MOIST. EXP. COEF. 

BTAHCl 

.4738-02 

27 


BTAHC2 

.1121+00 

28 


BTAHC3 

.1121+00 

29 

FLEXURAL MODULI 

EHCIF 

.1560+08 

30 


EHC2F 

. 1003+07 

31 

STRENGTHS 

SHC23 

.4601+04 

32 


SHCIF 

.1912+06 

33 


SHC2F 

.9411+04 

34 


SHCSB 

.6902+04 

35 

FIBER VOL. RATIO 

VFH 

.5500+00 




FLEXURAL STRENGTH 
***************** 


I })n|nK « 4c « }|c «« ^ )K :((«««« )|c «« infnic «««* 4c }|c 1c « )K «« )|u(c })nK }K * )K)K «««««« )(n|nK * 4nK ** 1 

I THEORY I LONGITUDINAL I TRANSVERSE I 

1 4c 4C 4C4C 4( 4c4( 4c 4c 4c 4c 4c 4c 4( 4( 4( 4( 4c 4t 4c4c 4c 4c 4( 4c 4c 4( 4( « 4( 4t 4c 4c 4c 4c 4c 4c 4c 4( 4( 4( 4c 4c 4c 4C 4t 4( 4c 4c4c 4c4c 4C4C4C 4( 4( 4c 4( 4( 4c 4cl 
I BILINEAR I .1530+06 I .7529+0^ I 

I SOLUTION I II 

1 4c 4c 4c 4c 4c 4c 4c 4C 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c I 
I LINEAR RECTANGULAR I .1882+06 I .8302+04 I 

I SOLUTION III 

1 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4C 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4C 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c I 
I NEUTRAL AXIS SHIFT I .1882+06 I .8802+04 I 

I SOLUTION I I I 

I4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4cX 
I PARABOLIC I .1912+06 I .9411+04 I 

I SOLUTION I I I 

1 4c 4c 4c 4 c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c I 
I RECTANGULAR I .2295+06 I .1129+05 I 

I SOLUTION III 

1 4c 4C 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4C 4C 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4C 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c I 


Using IATEST2 


This method is used to check the interactive mode of INHYD. IATEST2 is a dataset that, when 
called, automatically inputs data as the user would enter it point by point in the interactive mode. It 
can be used to test this mode quickly. This INHYD. IATEST2 dataset is shown here: 


0.20E+00 

PRIMARY 

FIBER 

AS GRAPHITE FIBER 

32.E+06»2.E+06f 2.E+06» 1 .E+06f .2f .25 
- .56E-06t 5.6E-06» .06» 1 . E+4 r 3 . E~04 f . 17 
580.0f 58.0r58 .Of 400.E+03f 400.E+03 
THAT 

EPOXY MATRIX 3501-5 - 70F - DRY 

.46E+06f . 1643E+06f .4f 32.E-06* .443E-01 r .25 
1 .25f 6.8E+03f 36.3E+03f7.E+03f 0. 4,2.0E~04 
273.0f 460.0f 294 .0 f 0 .00 
55/45 AS/E 
.55f . 45f 0.0 
HTM 

SECONDARY 

FIBER 

S-GLA3S FIBER 

12.4E+06f 12.4E+06f 5. 17E+06f 5. 17E+06f .2t .20 
2.8E-06f 2.8E“06f . 09 f 204 . , 3 . 6E-04 f . 17 
7.5f7.5f 7.5f .36E106f .36E+06 
THAT 

EPOXY MATRIX - 3501-5 •• 70F - DRY 

*46E+06f . 1643E+06f .4f 32.E-06f .443E-01 f .25 
1 .25f 6.8E + 03f 36.3E + 03f 7.E + 03f 0.4f 2.0E--04 
273.0t460.0f 294.0 f 0.00 
55/45 S-G/E 
♦55f .45f0.0 
HTM 

80/20 AS/E//S-G/E 
RULE OF MIXTURES 



To start a run, the user must follow this procedure, after logging on a terminal: 


@ASG. A INHYD 
©MAP, I .ABS 
IN INHYD. 

END 

©XQT .ABS 

The computer returns INHYD to be run in the batch or interactive (B/I) mode; the user selected 
interactive. Instead of entering the secondary composite volume ratio (which is in the dataset), the 
user enters ADD INHYD.IATEST2, as shown here. 

INHYD INTERACTIVE DATA ENTRY 

SECONDARY COMPOSITE VOLUME RATIO 
>eadd inhyd ♦ iatest2 


WHICH COMPOSITE DATA TO BE ENTERED <P/S) 
PROPERTIES TO BE ENTERED (C/F/M/T) 

INPUT DATA HEADING ( 1 TO 80 CHARACTERS) 
ENTER THE FIBER PROPERTY VALUES 


EFPl 

EFP2 

GFP12 

GFP23 

NUFP12 NUFP23 

CTEFPl 

CTEFP2 

RHOFP 

NFP 

DIFP CFPC 

KFPl 

KFP2 

KFP3 

SFPT 

SFPC 


PROPERTIES TO BE ENTERED (C/F/M/T) 


INPUT DATA HEADING ( 1 TO 80 CHARACTERS) 

ENTER THE MATRIX PROPERTY VALUES 

EMP GMP NUMP CTEMP RHOMP CMPC 

KMP SMPT SMPC SMPS BTAMP DIFMP 

ENTER THE THERMAL/MOISTURE PROPERTY VALUES 
TO TGDR TEMP MOIST (%) 

COMPOSITE HEADING (1 TO 80 CHARACTERS) 

ENTER CONSTITUENT VOLUME RATIOS 
FIBER MATRIX VOID 

METHOD FOR CALCULATING COMPOSITE PROPERTIES (H/M) 

WHICH COMPOSITE DATA TO BE ENTERED (P/S) 

PROPERTIES TO BE ENTERED (C/F/M/T) 

INPUT DATA HEADING (1 TO 80 CHARACTERS) 

ENTER THE FIBER PROPERTY VALUES 
EFSl EFS2 GFS12 GFS23 NUFS12 NUFS23 
CTEFSl CTEFS2 RHOFS NFS DIFS CFSC 

KFSl KFS2 KFS3 SFST SFSC 

PROPERTIES TO BE ENTERED (C/F/M/T) 

INPUT DATA HEADING ( 1 TO 80 CHARACTERS) 

ENTER THE MATRIX PROPERTY VALUES 

EMS GMS NUMS CTEMS RHOMS CMSC 

KMS SMST SMSC SMSS BTAMS DIFMS 

ENTER THE THERMAL/MOISTURE PROPERTY VALUES 
TO TGDR TEMP MOIST (Z) 
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COMPOSITE HEADING (1 TO 80 CHARACTERS) 
ENTER CONSTITUENT VOLUME RATIOS 



FIBER MATRIX VOID 


METHOD FOR CALCULATING COMPOSITE PROPERTIES (H/M> 

HYBRID HEADING <1 TO 80 CHARACTERS) 

METHOD FOR CALCULATING HYBRID PROPERTIES (R/M) 

The computer then tabulates the data taken from the dataset INHYD JATEST2. This tabulation is 
done in exactly the same format as is shown for the example with moisture and temperature effects 
used earlier and therefore is not listed here. The listing includes the tabulation of the hybrid 
composite properties and the flexural strengths. 


INHYD Batch Mode 

From Terminal 

The program is started as with the interactive mode according to the system in use. At the Lewis 
Research Center on the UNIVAC 1110, the following procedure is used: 

@ASG. A INHYD 
©MAP, I .ABS 
IN INHYD. 

END 

©XQT .ABS 
B 

The command selects the batch mode. Next the user adds a previously prepared dataset that contains 
the needed input data in the form desired. For example, 

©ADD INHYD.DATl 

The format of the input data for running INHYD in a batch mode is shown on the next page for 
INHYD.DATl. The “scale” entries were added to orient the input data in the proper columns and 
are not part of the actual data entry. Each table entry is identified after the input table. 


>0edfU inhyd.datl 
CASE UPPER ASSUMED 

ED 1 6R 1 A-FR I -04/22/83-13 : 59: 43- ( 23 » 24) 

EDIT 

0 ♦ >5cale 

123456789012345678901234567890123456789012345678901234567890123456789012 

o:>i 


0.20E+00 

i:>P 20 


0*20E+00 

PRIMARY 

FIBER 

AS GRAPHITE FIBER 
32*E+06 

2.E+06 

2.E+06 

1 .E+06 

^ 2 

*25 

-.56E-06 

5.6E-06 

.063 

1.EF04 

3.E-04 

.17 

580.0 

58.0 

58*0 

400*E+03 

TMAT 

EPOXY MATRIX - 3501 
.46E+06 

400.E+03 

-5 - 70F 

1643E+06 

- DRY 

* 4 

32.E-06 

.443E-01 

*25 

1*25 

14.0E+03 

36.3E+03 

16.0E+03 

0.4 

273*0 

2.0E-04 

460*0 

394*0 

0.00 


55/45 AS/E 

.55 

.45 

0*0 
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HTH 

SECONDARY 

FIBER 

S-GLASS FIBER 
20 ♦ >scale 

1234567890 1234567890 1234567890 12345678901 23456789012345678901 234567890 12 
20t>p 20 


S-GLASS FIBER 
12*4E+06 

12.4E+06 

5.17E+06 

5.17E+06 

.2 

%20 

2.8E-06 

2.8E-06 

.09 

204. 

3*6E-04 

.17 

7.5 

7.5 

7.5 

.36E+06 

THAT 

EPOXY MATRIX - 
♦46E+06 

.3E+06 

3501-5 - 70F 

.1643E+06 

- DRY 

.4 

32*E-06 

*443E-01 

%25 

1.25 

14.0E+03 

36.3E+03 

16.0E+03 

0.4 

273*0 

2.0E-04 

460.0 

394.0 

o 

o 

o 


55/45 S-G/E 

.55 

HTM 

♦ 45 

0.0 




80/20 AS/E//S-G/E 
ROM 

E0F:35 SCANtlS 
0 • >sc3le 


123456789012345678901234567890123456789012345678901234567890123456789012 

The entries in the DATl input table (symbols defined in appendix C) are as follows: 


PRIMARY 

FIBER 

AS GRAPHITE FIBER 


EFPl 

EFP2 

GFP12 

GFP23 

NUFP12 

NUFP23 

CTEFPl 

CTEFP2 

RHOFP 

NFP 

DIFP 

CFPC 

KFPl 

KFP2 

KFP3 

SFPT 

SFPC 




TMAT 





EPOXY MATRIX - 3501 

-5 - 70F - 

DRY 


EMP 

GMP 

NUMP 

CTEMP 

RHOMP 

CMPC 

KMP 

SMPT 

SMPC 

SMPS 

BTAMP 

DIFMP 




TO 

TGDR 

TEMP 

XMOISTURE 


55/45 A8/E 





FVR 

MVR 

VVR 



HTM 

SECONDARY 

FIBER 





B-GLASS FIBER 




EFSl 

EFS2 

GFS12 

GFS23 

NUFS12 

NUFS23 

CTEFSl 

CTEFS2 

RHOFS 

NFS 

DIFS 

CFSC 

KFSl 

KFS2 

KFS3 

SFST 

SFSC 




TMAT 





EPOXY MATRIX - 3501 

-5 - 70F - 

DRY 


EMS 

GMS 

HUMS 

CTEM3 

RHOMS 

CMSC 

KMS 

SMST 

SMSC 

SMSS 

BTAMS 

DIFMS 




TO 

TGDR 

TEMP 

•/MOISTURE 


55/45 S-G/E 





FOR 

MOR 

OVR 



HTM 

80/20 AS/E//S- 
ROM 

-G/E 





DATl contains the volume fraction of secondary composite as the first data entry; it also has the 
volumetric proportions of fiber, resin (matrix), and voids as well as all of the necessary fiber and 
matrix properties. All of these data are arranged in the proper format so that the computer can access 
the data and use it to calculate the primary and secondary composite properties. The program allows 
for temperature and moisture effects, so the subroutine HTM (which allows for hygrothermo- 
mechanical effects) is used to calculate the primary and secondary composite properties. The hybrid 
composite properties are then calculated by the rule of mixtures. Of course, any other property values 



(and headings) desired can be substituted for those presently contained in the sample data set to make 
it applicable to any desired intraply hybrid composite. 

DAT2 is to be used if primary composite properties are known but secondary composite properties 
must be calculated from fiber and matrix properties; otherwise, it is similar to DATl (see compiled 
listing). DAT3 uses fiber and matrix properties to calculate primary composite properties when 
secondary composite properties are known. DAT4 can be used when primary and secondary 
composite properties are known; it simply uses the rule of mixtures to determine hybrid composite 
properties. DATS is similar to DATl but with moisture present. Several other examples can be 
obtained by doing @ADD INHYD.DAT6, 7, 8, or 9 instead of @ADD INHYD.DATl as was done 
in the preceding example. Any other desired datasets can be made up by using the same format. 
Using INHYD. START results in a set of several sequential computer runs in which such things as 
primary and secondary composite ratios in the hybrid composite are varied or in which different 
properties for the matrix modulus (for example) are varied to study these effects on the calculated 
hybrid composite properties. To make INHYD. START run, follow this procedure: 

@ASG. A INHYD 
©MAP, I .ABS 

IN INHYD. 

END 

©ADD INHYD.START 


Initiated by Cards 

The cards needed to initiate an INHYD run at Lewis are 

“RUN, M/TP MSBJS, YOH1965, USERS NAME, 20, 200 

“Run card includes the project number, the time, and the number of pages sillowed but 
would vary with the requirements of the computer center. 

“COL 9000FD 
©ASG. A INHYD 
©MAP, I .ABS 
IN INHYD. 

END 

©XQT .ABS 
B 

©ADD INHYD.DATl 

It selects the specific run desired from those on file. 

©COL llOOFD 
“FIN 


Concluding Remarks 

This users manual provides step-by-step instructions for using the computer program INHYD 
(INtraply HYbrid composite Design). The instructions are given for both the user-interactive mode 
and the batch mode of running INHYD. The various options provided in INHYD are demonstrated 
with appropriate computer runs. Data input formats, sample computer outputs, and symbol lists are 
included. The general structure, capabilities, and options of INHYD are described in an appendix. 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio, August 12, 1983 



Appendix A 

Elements and Subroutines 

INHYD consists of a number of elements and subroutines; they are listed alphabetically and 

described briefly in this section 

ASD 17 properties of AS graphite in format to be used in INHYD 

CHANGE 1 Edits FIBMT and COMPP to remove comment symbols (C) from three cards 
(lines) to activate later commands; then calls INHYD.COM when FIBMT 
and COMPP are used 

CHANGE 2 Edits INHYD.FIBMT and COMPP, adding and removing C’s from several 
lines; then calls INHYD.COM 

CHANGE 3 Identical to CHANGE 1 

CHANGE 4 Adds and removes C’s before several lines of FIBMT and COMPP; then calls 
INHYD.COM 

CHANGE 5 Adds and removes C’s before several lines of FIBMT and COMPP; then calls 
INHYD.COM 

COM Packs and prepares INHYD; compiles FIBMT and COMPP; then initiates 

INHYD.START 

COMPILE Starting INHYD. COMPILE gives a compiled listing of the entire INHYD 

program (This depends on a UNIVAC control language program at Lewis 
and would not work at another computer center.) 

COMPP Generates hybrid composite properties based on micromechanics theory 

DATl Batch run of INHYD for properties of 80/20 AS/E//3501 at 250° F. Calcula- 

tions are made by using hygrothermomechanical (HTM) theory because of 
high-temperature use. 

DAT2 Calculates hybrid composite properties for 80/20 AS/E//S-Glass/E composite. 

The 37 primary composite properties for a 55/45 AS graphite/epoxy resin 
composite are read in. The secondary composite properties are calculated by 
FIBMT from the S-Glass fiber and 3501 epoxy matrix properties. Then the 
intraply hybrid composite properties are calculated by the rule of mixtures. 

DAT3 Same as DAT2 except that the secondary composite properties are read in and 

the primary composite is calculated by FIBMT from AS graphite and 3501 
resin properties. The rule of mixtures is used to calculate hybrid composite 
properties. 

DAT4 Calculates hybrid composite properties for 80/20 AS/E//S-Glass/E directly 

from the 37 primary composite and 37 secondary composite properties that 
are read in. The rule of mixtures is used. 

DAT5 Reads in fiber and matrix properties for both primary and second 2 iry 

composites. Program degrades matrix properties for the given moisture. 
Primary and secondary composite properties are calculated by using FIBMT 
(HTM can be substituted in place of FIBMT). Finally, the 80/20 AS/E//S- 
Glass/E hybrid composite properties are calculated from the primary and 
secondary composite properties by using COMPP, which uses the micro- 
mechanics theory. 

DAT6 Calculates the properties of a uniaxial fibrous composite from its fiber and 

matrix properties by using FIBMT 

DAT7 Calculates the hybrid composite properties for 80/20 HMS/E//S-Glass/E 

directly from the 37 primary composite and 37 secondary composite 
properties that are read in. The rule of mixtures is used. 
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DATS 


DAT9 


EPOXYD 

EPOXYMD 

EXE 

EXE2 

EXP 

FIBMT 

FLEXX 


HMSD 

HTM 

IATEST2 


IATEST3 

INSTR 

INTER 

KEVD 

LIST 


MAIN 


MIXRUL 

PRINT 


RUN 


Calculates the hybrid composite properties for 80/20 AS/E//S-Glass/E by 
using micromechanics theory. Input is the fiber and resin properties of the 
constituents. 

Calculates the hybrid composite properties for 80/20 AS/E//S-Glass/E. Input 
is the fiber and resin properties of the constituents. Primary and secondary 
composite properties are calculated by using micromechanics theory. Hybrid 
composite properties are then calculated by using the rule of mixtures. No 
temperatures or moisture effects are involved. 

Contains the 12 properties of an epoxy resin in proper format to be used by 
INHYD 

Same as EPOXYD but with a different value for resin modulus 

Contains commands in INHYD and END 

Identical to EXE 

Contains fiber and matrix properties for a 57/43 S-Glass/3501 composite. No 
moisture or temperature effects are involved. 

Calculates primary or secondary composite properties by using micro- 
mechanics theory 

Calculates and prints the flexural strength of any composite by using several 
theories (bilinear, linear rectangular, neutral axis shift, parabolic, and 
rectangular) 

Contains the 17 properties of HMS graphite fibers 

Calculates primary or secondary composite properties by using hygrothermo- 
mechanical theory 

Automatically inputs data as a user would enter it in interactive mode for 
80/20 AS/E//S-Glass/E hybrid composite; has temperature effect; moisture 
effect could be added. 

Automatically inputs data as a user would enter it in interactive mode for 
room temperature, dry, 80/20 AS/E//S-Glass/E hybrid composite 

Dataset containing instructions and explanations useful to INHYD user 

Interactively queries user for input data required for INHYD in interactive 
mode and then stores it for later use by main program 

Contains the 17 fiber properties of Kevlar 49 

To get a listing of INHYD, do ©START INHYD. LIST. This is not a 
compiled listing (see COMPILE). (This depends on a UNIVAC control 
language program at Lewis and it would not work at another computer 
center.) This command is used for listing and compiling convenience only. 

Computes intraply hybrid properties as well as primary and secondary 
composite properties. Output includes elastic moduli, shear moduli, Poisson’s 
ratios, thermal and moisture coefficients, and strengths. The equations are 
based on references 2 to 4 (appendix B). 

Predicts hybrid composite properties based on the rule of mixtures. See ref. 1 
(appendix B). 

This command prints the elements in the file INHYD.; also gives listings of 
the following INHYD subroutines: MAIN, HTM, FIBMT, FLEXX, 
MIXRUL, and COMPP. 

This command, used in the batch mode, starts a test run that includes 
INHYD. ST ART and INHYD. ST ART2 (see below); it also includes several 
datasets identified as CHANGES (1 to 5) (see above). Use care in invoking 
this command; it results in an output of about 3000 pages. 



RUNTESTl 

SGD 

START 

START2 

TAPE 


Starts INHYD.RUN from a terminal but gives listing from the computer high- 
speed printer (does not complete run) 

Has the 17 properties of S-Glass fibers needed for determining the properties 
of a S-Glass/resin composite 

Causes several sequential computer runs to be printed out of calculated prop- 
erties of hybrid composites with different primary and secondary composite 
ratios and composites with different fiber makeup 

Same as INHYD. START but uses a resin with different properties (modulus) 
than START 

List of INHYD elements in alphabetical order 
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COMPUTER CODE FOR INTRAPLY HYBRID COMPOSITE DESIGN* 


C. C. Chamis and J. H. Sinclair 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


ABSTRACT 

A computer program has been developed and is described herein for intra- 
ply hybrid composite design (INHYD). The program includes several composite 
micromechanics theories, intraply hybrid composite theories, and a hygrother- 
momechanical theory. These theories provide INHYD with considerable flexi- 
bility and capability that the user can exercise through several available 
options. Key features and capabilities of INHYD are illustrated through selec- 
ted samples. 


INTRODUCTION 

Intraply hybrid composites generally have two kinds of fibers embedded in 
the matrix. They have evolved as a structural material as a logical sequel to 
conventional composites and to interply hybrid composites. Intraply hybrid 
composites have unique features that can be used to meet diverse and competing 
design requirements in a more cost-effective way than either advanced or con- 
ventional composites. Some of the specific advantages of intraply hybrids 
over the constituent composites are balanced strength and stiffness, balanced 
bending and membrane mechanical properties, balanced thermal distortion sta- 
bility, reduced weight or cost, improved fatigue resistance, reduced notch 
sensitivity, improved fracture toughness or crack-arresting properties, and 
improved impact resistance. By using intraply hybrids, the designer can ob- 
tain a viable compromise between mechanical properties and cost to meet speci- 
fied design requirements for aerospace structures. 

The mechanical behavior of intraply hybrids has been investigated at the 
Lewis Research Center theoretically and experimentally (refs. 1 to 3). The 
theoretical methods and equations described in these references together with 
those for hygrothermal effects (ref. 4) have been integrated into a computer 
code for predicting hygral, thermal, and mechanical properties of, and thereby 
designing, intraply hybrid composites. This code is identified as INHYD for 
INtraply HYbrid-composite Design. The objective of this paper is to describe 
INHYD with respect to theory, equations, input, output, and its various 
options. 

The capability and various options of the program are briefly described. 
Typical input for INHYD includes fiber and resin properties or composite prop- 
erties, volume ratios of the primary and secondary fibers, glass transition 
temperature of the resin, moisture, and cure and use temperatures. Typical 
output includes moisture expansion coefficient, thermal expansion coeffi- 
cients, moduli and strengths (in plane and flexural). Selected samples of 


*This appendix is a reprint of NASA TM-82593. 



input and output are included to illustrate the flexibility and capability 
ofINHYD. Also, the planned extensions and couplings with integrated computer 
programs are briefly discussed. 


INTRAPLY HYBRID GEOMETRY AND DEFINITIONS 

The generic geometry of intraply hybrids included in INHYD is shown sche- 
matically in figure 1. The features to be noted in this figure are as follows: 

(1) The schematic shows an eight-ply, unidirectional intraply hybrid. 

Each ply consists of primary composite (blank) and secondary composite (cross 
hatched) . 

(2) The location of the secondary composite is regular and staggered 
through the thickness of the intraply hybrid. 

(3) The coordinate reference axes are 1, 2, and 3, where 1 is taken par- 
allel to the fiber direction, 2 is taken transverse to the fiber direction, 
and 3 is taken through the thickness. These coordinate axes will be called 
the material axes. All properties to be discussed in subsequent sections are 
defined with numerical subscripts corresponding to the material axes. 

It is clear from figure 1 that an intraply hybrid is made by laying tapes 
(tows) of primary composite and tapes (tows) of secondary composite in a defi- 
nite manner in order to obtain the desired volume ratios of each. For pur- 
poses of identification, the primary composite in the hybrid is the one that 
constitutes the largest volume ratio. This imposes no limitations on either 
the capability or flexibility of INHYD. Another way to view an intraply 
hybrid is the one in which the constituent composites are uniformly dispersed 
in a typical cross-sectional area. 


THEORIES INCLUDED IN INHYD 


The theories programmed in INHYD have been reported previously and in- 
clude (1) rule of mixtures (ref. 1), (2) micromechanics equations for hygro- 
thermal effects (ref. 2), and (3) intraply hybrid composite micromechanics 
(ref. 4). Also, several equations have been programmed for predicting flex- 
ural strength and through-the-thickness shear as will be described later. 

These theories consist of equations that predict intraply hybrid composite 
properties based on constituent composite (primary and secondary) properties. 
The constituent composite properties may in turn be predicted by micromechanics 
equations using fiber and matrix properties. Essential features of these 
theories are summarized below. The detailed derivations, equations, and jus- 
tifications are found in the original references. 

The rule-of-mixtures equations have the following form (ref. 1): 


*"hc ~ ''pc^pc '^sc’^sc 


( 1 ) 


where P denotes property and V volume ratio. The subscript HC denotes 
hybrid composite; PC, primary composite; and SC, secondary composite. Any 
property of the intraply hybrid composite can be predicted by using 
equation (1) when the properties of the primary and secondary composites have 
been specified. Equation (1), though of simple form, does satisfy the three 
principles of mechanics: force equilibrium, strain compatibility, and stress- 

strain relationships. Equation (1) predicts values that are in good agreement 
with measured data (ref. 1) and with other more sophisticated theories 
(ref. 2). Equation (1) is programmed as two subroutines in INHYD: one for 



the primary composite and one for the secondary composite. These subroutines 
require input of constituent composite properties. The individual properties 
required are summarized in later sections. 

The micromechanics equations from reference 4 are programmed in INHYD for 
predicting composite properties by using fiber and matrix properties and by 
accounting for hygrothermal effects. These equations are also programmed in 
two different subroutines: one for predicting primary composite properties 

and one for predicting secondary composite properties. Calls to these sub- 
routines require inputs of fiber and matrix properties, fiber volume ratio, 
void volume ratio, use temperature, and moisture. The individual properties 
needed are described later. 

The intraply hybrid composite micromechanics equations from reference 2 
are programmed in INHYD. Three subroutines are used: one each for the pri- 

mary and the secondary constituent composites and one that combines these two 
into an intraply hybrid. The first two subroutines require fiber and matrix 
properties; the third requires the primary and secondary composite properties 
predicted by the first two. 

The properties predicted by all of these theories are those that result 
from uniform response (due to moisture, temperature, or stress) variation 
through the thickness of the intraply hybrid. Properties resulting from non- 
uniform response variation through thickness such as flexural strength and 
short-beam shear strength are considered in the next section. 


FLEXURAL PROPERTIES 

The equations programmed for flexural properties in INHYD are for flex- 
ural strength along the 1- and 2-directions and the accompanying shear strength 
in the 1-3 plane (short beam shear) in figure 1. Four different equations are 
programmed for flexural strengths. The equations have not been reported pre- 
viously. The equations are readily derivable using the well-known simple beam 
theory and assumed stress variations through the beam thickness at fracture 
(fig. 3). The four equations included in INHYD were derived by assuming 
(1) bilinear stress distribution; (2) parabolic stress distribution; (3) 
linear rectangular stress variation; and (4) rectangular stress distribution 
(fig. 3). All these distributions are possible in intraply hybrid composites 
depending on their constituents. The selections of a specific equation re- 
quires user interaction based on measured data correlation. The default op- 
tion in INHYD is the parabolic stress distribution. All four equations are 
programmed in two different subroutines: one for the primary composite and 

one for the secondary composite. The corresponding flexural strengths for tne 
intraply hybrid are calculated in the same subroutine for intraply hybrid 
micromechanics. 

The equation for the through-the-thickness shear strength in the 1-3 plane 
is derived by assuming parabolic shear stress variation. The integral of this 
parabolic variation is equated to uniform in-plane shear to obtain the desired 
equation. This results in the following simple relationship: The through-the- 

thickness shear (short-beam shear) is equal to 1.5 times the corresponding 
in-plane shear strength. The short-beam shear strength of the intraply hybrid 
is calculated in the same subroutine as the flexural strength. The equation 
for predicting the short-beam shear strength in the 2-3 plane (fig. 1) is sim- 
ilar, but it is not programmed in INHYD. The reason that it is not programmed 
is that experimental data indicate that specimens tested for this shear fail 
invariably by transverse flexure. 



INHYD COMPUTER PROGRAM STRUCTURE 


The logic flow of the INHYD computer program is illustrated schematically 
in figure 2. The theoretical functions of the program are enclosed with double 
lines. The input, intermediate output, and final output are enclosed with 
single lines. Typical final uses of the output are enclosed with dashed 
lines. The various micromechanics theories and requisite subroutine described 
previously are included in the first double-line block. Those for predicting 
the intraply hybrid properties are included in the second double-line block. 

The user controls INHYD through the main subroutine by selecting combinations 
of the several available options. The options in INHYD are summarized in 
table I. It is clear from table I that combinations of options can be speci- 
fied for theory, hybrid, input property sets, and hygrothermal effects. 

The input of INHYD is summarized in table II. A sample user input data 
set when the intraply hybrid properties will be predicted by using fiber and 
matrix properties is shown in table III . The first set of Booleans are the 
options for input read-in formats and for environmental effects. The first 
line of data is the secondary composite volume ratio. The second set of 
Booleans specify the composite micromechanics subroutines to be used. The 
alphanumerics describe the primary composite system. The next (third) set of 
Booleans define the fiber properties that will be read in. This set of 
Booleans is followed by the primary composite fiber properties: moduli, 

Poisson's ratios, thermal expansion coefficients, density, number of fibers 
per end, fiber diameter, heat capacity, heat conductivities, and tensile and 
compressive strengths. The alphanumeric card following the fiber properties 
describes the primary composite matrix and the environmental conditions. The 
next (fourth) set of Booleans specify the matrix properties to be read in. 

These matrix properties are moduli, moisture expansion coefficient per percent 
of moisture, thermal expansion coefficients, heat capacity, thermal conductiv- 
ity, strengths (tensile, compressive, and shear), moisture diffusivities, glass 
transition temperatures, reference temperature, moisture content, and volume 
ratios (fiber, matrix, and void). The data following are the input data needed 
for the fiber and matrix in the secondary composite. The format for these data 
is similar to those (one for one correspondence) for the primary composite. 

Formatted output (compiled data) of the input data for the primary com- 
posite fiber properties is shown in table IV and in table V for the matrix. 

In these tables the Booleans for the various program options, the fiber and 
matrix types, the volume ratios, and the constituent properties (with the cor- 
responding program name) are displayed. The reader can readily match corres- 
ponding values from the user input data (table III) and the compiled data 
(tables IV and V). 

INHYD program output features are summarized in table VI. A sample out- 
put for the predicted primary composite properties (using the input data 
tables IV and V) is illustrated in table VII. The volume ratio values and 37 
properties are listed in this table. The various properties are identified by 
both name and assigned variable used in the program. Note the flexural prop- 
erties, lines 29 to 34. Note also the ply thickness (line 35), the interply 
thickness (line 36), and the interfiber spacing (line 37). The format for the 
outputs of the secondary composite and for the intraply hybrid are similar to 
that for the primary (table VII). 
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INHYD PROGRAM CAPABILITY - ILLUSTRATED EXAMPLES 


The INHYD program capability and flexibility is illustrated, in part, 
with typical examples. The effects ^of hybridization on selected properties 
are summarized in table VIII for 70° F, dry environmental conditions. The 
properties shown in this table were predicted by using the user input data in 
table III. Corresponding properties for the primary, the secondary, and the 
intraply hybrid composite are summarized. Corresponding properties for the 
primary composite are the same as those in table VII. 

The thermal degradation, due to high temperature, in the properties of 
two intraply hybrids (AS/E//S-G/E and (AS/E//KEV/E) is illustrated in table 
IX. It is instructive to note in this table that the thermal degradation is 
insignificant in longitudinal properties (except longitudinal compressive 
strength). However, the thermal degradation is considerable (about 30 percent) 
for transverse and shear properties and for longitudinal compressive strength. 

The hygrothermal degradation in the properties of the AS/E//S-G/E intra- 
ply hybrid is illustrated in table X. The hygrothermal environments are 70° F 
with 1 percent moisture and 250° F with 1 percent moisture. Again, the hygro- 
thermal degradation is insignificant for the longitudinal properties, except 
for the longitudinal compressive strength. On the other hand the 1 percent 
moisture degrades the transverse and shear properties and the longitudinal 
compressive strength about 10 percent for the 70° F case and about 30 percent 
for the 250° F case. 


INHYD EXTENSIONS AND COUPLING 

The properties predicted by INHYD described previously constitute only the 
first part of INHYD. The program can be readily extended to predict other 
properties of intraply hybrids such as impact resistance and fatigue. In addi- 
tion it can be made a part of or coupled with integrated programs for special 
or general structural analysis. The near-future planned extensions and cou- 
plings with other programs are summarized in table XI. The theory for impact 
resistance in reference 5 will be used. The theory for fatigue resistance is 
presently under development and will be reported in the ASTM 6th Conference on 
Composite Materials: Testing and Design, May 1981. INHYD will be coupled with 

an in-house laminate analysis code (MFCA, ref. 6). It will also be coupled 
with three integrated computer programs under in-house development: CODSTRAN - 

Composite Durability Structural Analysis (ref. 7); COBSTRAN - Composite Blade 
Structural Analysis (ref. 8); and CISTRAN - Composite Impact Structural 
Analysis (ref. 9). 


CONCLUSIONS 

A computer program called INHYD has been developed to predict the prop- 
erties of unidirectional intraply hybrid composites and, therefore, to assist 
in the design of these hybrids. Several composite micromechanics and intraply 
hybrid theories and a hygrothermal mechanical theory in INHYD provide the pro- 
gram with considerable flexibility that the user exercises through combina- 
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tions of options. These options control the input data, fiber and matrix, 
unidirectional composite or combinations, the output, and the theory. Selected 
samples illustrate key features and capabilities of INHYD. INHYD provides the 
designer or analyst with a convenient analytical means to investigate several 
intraply hybrids during the preliminary design phases. INHYD can also be used 
to guide, and therefore keep to a minimum, required characterization of intra- 
ply hybrids. 
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Figure 2. - Flow chart of INHYD computer program. 
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Figure 3. - Geometry, force, and stress distributions for estimating flexural strength. 


TABLE I. - PROGRAM OPTIONS 
Selects theory to be used 
Hybrid or conventional composite 
Types of property inputs 
Composite properties 
Constituent properties 
Combinations 
Hygrothermal effects 


TABLE II. - PROGRAM INPUT 

User options 

Constituent properties (fiber, matrix) or properties 
of primary and secondary composites or combinations 
Fiber volume ratio, void volume ratio 
Primary composite ratio, secondary composite ratio 
Temperature, moisture 
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TABLE III. - USER INPUT DATA 


T T T 

0.20E+00 

T 

T 

T 

T 

AS GRAPHITE FIBER 
F 
T 
F 
F 


32.EF06 
.25 
3.E-04 
400.EF03 
EPOXY HATRIX - 
F 
F 
F 
T 

.A6t^06 

.25 

0.4 

460.0 

.55 

F 

S-BLASS FIBER 
F 
T 
F 
F 


2.E+06 
~.56E~06 
.17 
400.Ef03 
3501-5 - 250F 


.1643E+06 

1.25 

2.0E-04 

394.0 

.45 


12.4E+06 

.20 


12.4E+06 

2.8E-06 


3.6E-04 .17 

.36E+06 .3E+06 

EPOXY MATRIX - 3501-5 - 250F 


2.E+06 

5.6E-06 

580.0 


DRY 


♦ 4 

6.8E+03 


0.0 

0.0 


5. 17E+06 
2.8E-06 
7.5 


l.E+06 

.063 

58.0 


32.E-06 

36.3E+03 


5.17E+06 

.09 

7.5 


l.E+04 

58.0 


.443E-01 

7.E+03 


2oit 

7.5 


DRY 



TABLE IV. - COMPILED DATA - FIBER 


vsc= 

« 200+00 





MRUL = 

T 

PHTMP= 

T 

SHTHF‘= 

T 

PRIM = 

T 





COMP = 

F 

FIBER= 

T 

HAT = 

F THAT* 


PRIMARY FIBER PROPERTIES 
AS GRAPHITE FIBER - 250F “ 


1 

ELASTIC MODULI 

EFPl 

.3200+08 

2 


EFP2 

.2000+07 

3 

SHEAR MODULI 

6FP12 

.2000+07 

4 


GFP23 

. 1000+07 

5 

POISSON'S RATIO 

NUFP12 

.2000+00 

6 


NUFP23 

.2500+00 

7 

THERM. EXP. COEF. 

CTEFPl 

.5600-06 

8 


CTEFP2 

5600-05 

9 

DENSITY 

RHOFP 

.6300-01 

10 

NO. OF FIBERS/END 

NFP 

. 1000+05 

11 

FIBER DIAMETER 

DIFP 

.3000-03 

12 

HEAT CAPACITY 

CFPC 

. 1700+00 

13 

HEAT CONDUCTIVITY 

KFPl 

.5800+03 

14 


KFP2 

.5800+02 

15 


KFP3 

.5800+02 

16 

STRENGTHS 

SFPT 

.4000+06 

17 


SFPC 

.4000+06 


TABLE V. - COMPILED DATA - MATRIX 


COMP= F FIBER= 

F 

MAT= F 

TMAT= T 


FRIHARY HATRIX PROPERTIES 
EPOXY MATRIX - 3S01-5 - 250F 

- DRY 




TGDR- .4600'f03 

T = 

.3940+03 

M= 

.0000 


ORIGINAL MATRIX PROPERTIES 


1 

ELASTIC MODULUS 

EMP 

.4600+06 

2 

SHEAR MODULUS 

GMP 

. 1643+06 

3 

POISSON'S RATIO 

NUMP 

.4000+00 

4 

THERM. EXP. COEF. 

CTEMP 

.3200-04 

5 

DENSITY 

RHOMP 

.4430-01 

6 

HEAT CAPACITY 

CMPC 

.2500+00 

7 

HEAT CONDUCTIVITY 

KMP 

.1250+01 

8 

STRENGTHS 

SMPT 

.6800+04 

9 


SMPC 

.3630+05 

10 


SMPS 

.7000+04 

11 

MOISTURE COEF. 

STAMP 

.4000+00 

12 

DIFFUSIVITY 

DIFMP 

.200C-03 
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TABLE VI. - OUTPUT 


• PRIMARY AND SECONDARY COMPOSITE PROPERTIES 

3 EACH, ELASTIC MODULI, SHEAR MODULI, POISSON'S RATIOS, THERM. EXP. COEF. 
DENSITY, HEAT CAPACITY 

3 HEAT CONDUCTIVITIES, 5 IN-PLANE STRENGTHS 

3 EACH, MOISTURE DIFFUSIVITIES AND MOISTURE EXPANSION COEFS. 

2 FLEXURAL MODULI 

4 THROUGH-THE-THICKNESS STRENGTHS 

PLY THICKNESS, INTERPLY THICKNESS, INTERFIBER SPACING 

• INTRAPLY HYBRID COMPOSITE PROPERTIES 

SAME AS ABOVE EXCEPT FOR LAST LINE PLUS COMPOSITE FIBER VOLUME RATIO 


TABLE VII. - TYPICAL COMPUTER OUTPUT OF INTRAPLY HYBRID PROPERTIES 

PRIHAKY COMPOSITE PROPERTIES 


• 

5500+00 VMP= 

.4500+00 

VVP» .0000 

1 

ELASTIC HODULI 

EPCl 

.1780+08 

2 


EPC2 

. 1034 + 07 

3 


EPC3 

.1034+07 

4 

SHEAR hODULI 

GPC12 

.4902+06 

5 


GPC23 

.4148+06 

6 


GPC13 

.4902+06 

7 

POISSON'S RATIO 

NUPC12 

.2900+00 

8 


NUPC23 

.2465+00 

9 


NUPC13 

.2900+00 

10 

THERM. EXP. COEF. 

CTEPCl 

1816-06 

11 


CTEPC2 

.1640-04 

12 


CTEPC3 

. 1640-04 

13 

DENSITY 

RHOPC 

.5458-01 

14 

HEAT CAPACITY 

CPC 

.2048+00 

15 

HEAT CONDUCTIVITY 

KPCl 

.3196+03 

16 


KPC2 

.3916+01 

17 


KPC3 

.3916+01 

18 

STRENGTHS 

SPCIT 

.2224+06 

19 


SPCIC 

. 1712 + 06 

20 


SPC2T 

.6407+04 

21 


SPC2C 

.3420+05 

22 


SPC12 

.6595+04 

23 

MOIST. DIFFUSIVITY 

DPCl 

.9000-04 

24 


DPC2 

.5168-04 

25 


DPC3 

.5168-04 

26 

MOIST. EXP. COEF. 

BTAPCl 

.4384-02 

27 


DTAPC2 

•1131+00 

28 


BTAPC3 

. 1131 + 00 

29 

FLEXURAL MODULI 

EPCIF 

. 1780 + 08 

30 


EPC2F 

.1034+07 

31 

STRENGTHS 

SPC23 

.6595+04 

32 


SPCIF 

.2418+06 

33 


SPC2F 

.1349+05 

34 


SPCSB 

.9893+04 

35 

PLY THICKNESS 

TPC 

.5000-02 

36 

INTERPLY THICKNESS 

PLPC 

.5850-04 

37 

INTERFIDER SPACING 

PLPCS 

.5850-04 
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TABLE VIII. - EFFECT OF HYBRIDIZATION ON COMPOSITE PROPERTIES 


[Dry; temperature, 70° F.] 


Property 

AS/E primary 
composite, 
80 percent 

S-G/E secondary 
composite, 

20 percent 

Intraply hybrid 
composite, 80/20 
AS/E//S-G/E 

Elastic modul i : 
EHCl 
EHC2 

0.1780+08 

0.1034+07 

0.7015+07 

0.1524+07 

0.1564+08 

0.1132+07 

Shear moduli: 
GHC12 
GHC13 

0.4902+06 

0.4148+06 

0.5517+06 

0.5517+06 

0.5025+06 

0.4422+06 

Poisson's ratio, NUHC12 

0.2900+00 

0.2900+00 

0.2900+00 

Strengths: 

SHCIT 

SHClC 

SHC2T 

SHC2C 

SHC12 

0.2224+06 

0.1712+06 

0.6407+04 

0.3420+05 

0.6595+04 

0.2037+06 

0.1448+06 

0.6407+04 

0.3420+05 

0.6595+04 

0.2187+06 

0.1659+06 

0.6407+04 

0.3420+05 

0.6595+04 

Flexural moduli: 
EHCIF 
EHC2F 

0.1780+08 

0.1034+07 

0.7015+07 

0.1524+07 

0.1564+08 

0.1132+07 


TABLE IX. - ILLUSTRATION OF THERMAL DEGRADATION IN INTRAPLY HYBRIDS 


Intraply hybrid com- 
posite property 

80/20 AS/E//S-G/E 

80/20 AS/E//KEV/E 

Temperature, “F 


70 

250 

70 

250 

Elastic modul i : 
EHCl 
EHC2 

0.1564+08 

0.1132+07 

0.1557+08 

0.8067+06 

0.1670+08 

0.9365+06 

0.1662+08 

0.6995+06 

Shear moduli: 
GHC12 
GHC23 

0.5025+06 

0.4422+06 

0.3368+06 

0.3077+06 

0.4504+06 

0.3715+06 

0.3104+06 

0.2693+06 

Poisson's ratio, NUHC12 

0.2900+00 

0.2900+00 

0.3065+00 

0.3065+00 

Strengths: 

SHCIT 

SHClC 

SHC2T 

SHC2C 

SHC12 

0.2187+06 

0.1659+06 

0.6407+04 

0.3420+05 

0.6595+04 

0.2175+06 

0.1432+06 

0.4040+04 

0.2157+05 

0.4159+04 

0.2227+06 

0.1531+06 

0.6407+04 

0.3420+05 

0.6595+04 

0.2217+06 

0.1306+06 

0.4040+04 

0.2157+05 

0.4159+04 

Flexural moduli: 
EHCIF 
EHC2F 

0.1564+08 

0.1132+07 

0.1557+08 

0.8067+06 

0.1670+08 

0.9365+06 

0.1662+08 

0.6995+06 




TABLE X. - ILLUSTRATION OF HYGROTHERMAL DEGRADATION IN 


INTRAPLY HYBRID (80/20 AS/E//S-G/E) 


Intraply hybrid com- 
posite property 

Temperature, “F 

70 

250 


Dry 

1,0 Percent 
moisture 

Dry 

1,0 Percent 
moisture 

Elastic moduli: 
EHCl 
EHC2 

0.1564+08 

0.1132+07 

0.1562+08 

0.1033+07 

0.1557+08 

0.8067+06 

0.1552+08 

0.5713+06 

Shear moduli : 
GHC12 
GHC23 

0.5025+06 

0.4422+06 

0.4500+06 

0.4006+06 

0.3368+06 

0.3077+06 

0.2279+06 

0.2140+06 

Poisson's ratio, NUHC12 

0.2900+00 

0,2900+00 

0.2900+00 

0.2900+00 

Strengths: 

SHCIT 

SHClC 

SHC2T 

SHC2C 

SHC12 

0.2187+06 

0.1659+06 

0.6407+04 

0.3420+05 

0.6595+04 

0.2183+06 

0.1584+06 

0.5626+04 

0.3003+05 

0.5791+04 

0.2175+06 

0.1432+06 

0.4040+04 

0.2157+05 

0.4159+04 

0.2169+06 

0.1297+06 

0.2631+04 

0.1405+05 

0.2709+04 

Flexural moduli: 
EHCIF 
EHC2F 

0.1564+08 

0.1132+07 

0.1562+08 

0.1033+07 

0.1557+08 

0.8067+06 

0.1552+08 

0.5713+06 


TABLE XI. - PLANNED EXTENSION 

• PROGRAM AND DOCUMENTATION HILL BE MADE AVAILABLE TO THE PUBLIC THROUGH 
COSMIC 

• EXTEND PROGRAM TO DETERMINE IMPACT RESISTANCE, DEFECT PROPAGATION, AND 
FATIGUE RESISTANCE OF UNIDIRECTIONAL INTRAPLY HYBRIDS 

• COUPLE WITH MULTILAYERED FILAMENTARY LAMINATE ANALYSIS (MFCA) 

• COUPLE INHYD WITH 

• CODSTRAN - COMPOSITE DURABILITY STRUCTURAL ANALYSIS 

• COBSTRAN - COMPOSITE BLADE STRUCTURAL ANALYSIS 

• CISTRAN - COMPOSITE IMPACT STRUCTURAL ANALYSIS 

• CONVERT PROGRAM TO IBM 370/3033 




Appendix C 
Symbols 

This appendix defines the symbols used in the computer program INHYD. Part I displays and 
defines the property symbols as they appear in INHYD. Part II lists in alphabetical order additional 
symbols that are used to represent subscripts. In the program the symbols from parts I and II are used 
together to further define the property. 


INHYD. Property Symbols 

Primary fiber properties: 

EFP1,EFP2 

GFP12,GFP23 

NUFP12,NUFP23 

CTEFP1,CTEFP2 

RHOFP 

NFP 

DIFP 

CFPC 

KFP1,KFP2,KFP3 

SFPT,SFPC 


elastic moduli 

shear moduli 

Poisson’s ratios 

thermal expansion coefficients 

density 

number of fibers per end 
fiber diameter 
heat capacity 
heat conductivities 
strengths 


Primary matrix properties: 

EMP 

GMP 

NUMP 

CTEMP 

RHOMP 

CMPC 

KMP 

SMPT,SMPC,SMPS 

STAMP 

DIFMP 


elastic modulus 

shear modulus 

Poisson’s ratio 

thermal expansion coefficient 

density 

heat capacity 

heat conductivity 

strengths 

moisture coefficient 
diffusivity 


Primary composite properties: 


EPC1,EPC2,EPC3 

GPC12,GPC23,GPC13 

NUPC12,NUPC23,NUPC13 

CTEPCl ,CTEPC2,CTEPC3 

RHOPC 

CPC 

KPC1,KPC2,KPC3 

SPC1T,SPC1C, 

SPC2T,SPC2C, 

SPC12 

DPC1,DPC2,DPC3 
BTAPCl ,BTAPC2,BTAPC3 


elastic moduli 

shear moduli 

Poisson’s ratios 

thermal expansion coefficients 

density 

heat capacity 

heat conductivities 

strengths 


moisture diffusivities 
moisture expansion coefficients 


EPC1F,EPC2F 

SPC23,SPC1F 

SPC2F,SPCSB 

TPC 

PLPC 

PLPCS 

Secondary fiber properties: 

EFS1,EFS2 

GFS12,GFS23 

NUFS12,NUFS23 

CTEFS1,CTEFS2 

RHOFS 

NFS 

DIFS 

CFSC 

KFS1,KFS2,KFS3 

SFST,SFSC 

Secondary matrix properties: 

EMS 

GMS 

NUMS 

CTEMS 

RHOMS 

CMSC 

KMS 

SMST,SMSC,SMSS 

BTAMS 

DIFMS 

Secondary composite properties, 

ESC1,ESC2,ESC3 
GSC12,GSC23,GSC13 
NUSC12,NUSC23,NUSC13 
CTESC 1 ,CTESC2,CTESC3 
RHOSC 
CSC 

KSC1,KSC2,KSC3 

SSC1T,SSC1C, 

SSC2T,SSC2C, 

SSC12 

DSC1,DSC2,DSC3 
BTASC1,BTASC2,BTASC3 
ESC1F,ESC2F 
SSC23,SSC1F, 


flexural moduli 
strengths 

ply thickness 
interply thickness 
interfiber spacing 

elastic moduli 

shear moduli 

Poisson’s ratios 

thermal expansion coefficients 

density 

number of fibers per end 
fiber diameter 
heat capacity 
heat conductivities 
strengths 

elastic modulus 

shear modulus 

Poisson’s ratio 

thermal expansion coefficient 

density 

heat capacity 

heat conductivity 

strengths 

moisture coefficient 
diffusivity 

elastic moduli 
shear moduli 
Poisson’s ratio 

thermal expansion coefficients 

density 

heat capacity 

heat conductivities 

strengths 

moisture diffusivities 
moisture expansion coefficients 
flexural moduli 
strengths 



SSC2F,SSCSB 

TSC 

PLSC 

PLSCS 

Hybrid composite properties: 

EHC1,EHC2,EHC3 

GHC12,GHC23,GHC13 

NUHC12,NUHC23,NUHC13 

CTEHCl ,CTEHC2,CTEHC3 

RHOHC 

CHC 

KHC1,KHC2,KHC3 

SHC1T,SHC1C, 

SHC2T,SHC2C, 

SHC12 

DHC1,DHC2,DHC3 
BTAHCl ,BTAHC2,BTAHC3 
EHC1F,EHC2F 

SHC23,SHC1F, 

SHC2F,SHCSB 

VFH 


ply thickness 
interply thickness 
interfiber spacing 

elastic moduli 
shear moduli 
Poisson’s ratio 

thermal expansion coefficients 

density 

heat capacity 

heat conductivities 

strengths 

moisture diffusivities 
moisture expansion coefficients 
flexural moduli 
strengths 

fiber volume ratio 


Primary matrix properties with 
TO 

TEMP 

KMST 

RHOMST 

TGDR 

MOIST<^o 

CMST 


moisture: 

reference temperature 

test temperature 

moisture conductivity 

moisture density 

dry glass transition temperature 

moisture content, percent 

moisture heat capacity 


INHYD. Subscripts 

B 

C 

F 

H 

P 

S 

1 

2 

3 


bending 

composite; compression 
fiber; flexural 
hybrid 
primary 

secondary; shear 

material axes parallel to fiber direction 
material axes transverse to fiber direction 
material axes through thickness 
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